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We study a one-dimensional elliptic problem with highly oscillatory random diffu- 
PLh ' sion coefficient. We derive a homogenized solution and a so-called Gaussian corrector. 

. We also prove a "pointwise" large deviation principle (LDP) for the full solution and 

(-H I approximate this LDP with a more tractable form. Applications to uncertainty quan- 

■ tification are considered. 
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1 Introduction 

Partial differential equations whose (deterministic or random) coefficients have fine- 
scale structure are notoriously difficult to solve. Here we consider the following elliptic 
problem: Given probability space {Q.,S,P), find Ue G -L^(i7; ifg ([0, 1])) satisfying 

—^A{x,^,^,uj)^Ue = f{x), ajG(0,l), WGO, £<1. (1) 

Here ^ : O — > is a random vector. In our case we assume either 

(i) A,{x) = A{x,^, f ,a;) = a{x,0 + bi^,u;) 

(ii) As{x) = ^(f ), with A{y)~^ = ^^y^ for a stationary dependent sequence 7^ (which 
depends on ^). [-J denotes the "floor" function. 

Case (i) is motivated by a Karhuncn-Locve expansion, which for Q : [0, 1] x — >■ M 
having continuous covariance T : [0, 1]^ M, T{x,y) := KQ{x)Q{y), is given by 
00 „ 
Q{x,oj) = A/Aj'0j(a;)/tj(x), where E{i/;iipj}= / hi{x)hj{x) dx = Sij. 
j=i 

See [18, 16]. This is multi-scale, but could be approximated by two scales. In case 
(i), the high frequency randomness f is "decoupled from in the sense that, after 
conditioning on ^, A{x,£,,y,uj) is stationary in y. So in the model ^ = a + 6 we 
are assuming 6 is a stationary random field. For simplicity we always assume f{x) is 
deterministic. 

Case (ii) is an example of dependent media with "short range correlations." Use of 
different kernels hk allows some flexibility in modeling. Although we only consider the 
case where the final jj are stationary, generalizations (using e.g. hj{x)) would not be 
difficult. 

In some PDE of interest (e.g. (1), or other elliptic equations [1], and also linear 
transport [4]) the solution admits an expansion of the form 

Ue{x) = Uq{x) + Ve{x) + Re{x), (2) 

where 77,0 is a homogenized solution that is dominant in the limit e — t- [17, 20], the 
remainder is negligible in some sense, and is given by an oscillatory integral 

v,{x) = J G{x,y)q{y,^)dy. (3) 
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In these cases one expects (and can often prove) that — uo) converges in dis- 

tribution to a Gaussian process (.x) (often a = d/2). Thus, one is justified in 
approximating by uq plus a Gaussian corrector: 

Uo + e°'v{x). (4) 

Prom an uncertainty quantification (UQ) perspective, this represents a significant sim- 
phfication. Computation of the homogenized solution uq is much less expensive than 
Ug. The corrector v{x) has an explicit form in terms of e.g. an Ito integral. This allows 
explicit calculation of the correlation function. Moreover, draws from the random pro- 
cess £°'v{x) can be done with minimal (compared to calculation of u^) effort. Another 
utility of corrector results is for validation of numerical homogenization schemes. For 
example, it is known that the numerical homogenization techniques MsFEM and HMM 
give solutions that converge to the correct homogenization limit uq as e, — t- 0. The 
question as to whether £~"(m^ — uq) converges to the correct limit is explored (for (1)) 
in [3]. 

As a downside, central limit approximations such as (4) arc expected to work well 
only for moderate deviations, i.e. for {ug — uo\ ^ 0(e°). Sometimes of interest in UQ 
applications are questions related to large deviations, e.g. P[u£ > £] for some £ ^ 0(1). 

Our main contribution is an investigation of the large-deviation behavior of Ug, 
the solution to (1). As it turns out, it is possible to derive a large deviation prin- 
ciple (LDP) for the solution Ue{x) (for fixed x). This gives asymptotic limits of e.g. 
elog P{u£{x) > £] for £ ~ 0(1). The resultant rate function is given implicitly (see the- 
orems 4.4, 4.5) as a solution to two (one convex and one non-convex) four-dimensional 
optimization problems. We also derive an approximate LDP (proposition 5.1) that cor- 
responds to the approximation u^ ~ uq+Vs- Since Vs is given explicitly as an oscillatory 
integral, the rate function is "more explicit" , being the result of a one-dimensional con- 
vex optimization problem. We verify numerically that the approximate rate function 
works well when £ <^ 1 and £ ~ 0(1) but "not too large." This sort of approximate 
LDP should be available in other situations where the solution can be approximated 
by a homogenized term and an oscillatory integral. Along the way we also derive a 
large deviation principle for some one dimensional oscillatory integrals, which appears 
to be new as well. 

A secondary contribution is a generalization of homogenization and corrector re- 
sults. Homogenization results typically start with a uniformly (over all realizations) 
elliptic diffusion coefficient of the form A{^,uj). In this case the homogenized tensor 
is constant. Here we generalize these results slightly by allowing for non-constant low- 
frequency randomness (in the case (i)) and relaxing the uniform ellipticity requirement 
to (8). We do this by conditioning on the coarse-scale and bounding higher moments 
of Ag'^. These assumptions are more in line with those encountered in UQ. We also 
prove almost sure convergence of the homogenized tensor. This is motivated by the fact 
that in practice the homogenized tensor can be obtained by picking one high-frequency 
media realization and averaging over a domain of size p [8]. Thus, it is nice to know 
that with probability one this realization converges as p ^ oo. 

Our large deviation result allows us to determine (roughly) to what degree the 
Gaussian corrector captures the tail behavior of the solution. This is useful in applica- 
tions where one may consider replacing the full solution with the homogenized solution 
plus a Gaussian corrector. Also, although we don't answer this here, questions such 
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as "does HMM capture the large-deviation behavior of the solution" could potentially 
be answered in a manner similar to the question "does HMM capture the moderate 
deviation behavior" as discussed above. Also of interest (and also not pursued further 
here) is the relation between large deviations and importance sampling of rare events. 
As it turns out, a large deviations result can give "asymptotically efficient" importance 
functions [10, 14]. 

In section 2 an asymptotic expansion of is presented. In section 3, theorems 
3.1 and 3.2 give results on the homogenized convergence Ue uo, and the corrector 
characterizing moderate deviations of — uq. In section 4 large deviations are con- 
sidered. After introducing the subject we derive a large deviation principle for two 
types of media, each with piecewise constant high-frequency parts (theorems 4.4, 4.5). 
We next present our approximate LDP in 5.1 and then numerical results in section 
6. Proofs of the homogenization and corrector theorems, which are generalizations of 
known results, are relegated to sections 7.1 and 7.2. 



2 Asymptotic expansion of the solution Us 

The boundary value problem (1) may be integrated and the solution is 

nix)-- r^d.+i^ r^ds (5) 

Here we define 1-average (•) of a function (f) : [0, 1] M, and F : [0, 1] ^ M by 



f mdy, F{s):= r f{t)dt. 
Jo Jo 



Define the homogenized tensor Aq by 

Ao{x): = {E^{A{x,0)-^})-\ 
where for X : ^ R, we denote conditional expectation by 

E^{X}: = E{X|e}= / X{u)dF^{u). 

Above the measure P^ = -P[-| is defined implicitly. 
Now re-write the integrals appearing in (5) as 



Ae{s) Jo Ao{s) ^' Jo Ae{s) Jq Ao{s) 

Defining the homogenized solution uq by the equation (5) with Aq rather than A^ (or 
equivalently the weak solution to (1) with Aq rather than A^), we have the following 
expansion. 

Us{x) = Uq{x) + Ve{x) + Re{x), 



v,{x) . - Y^ + [Y, X, ^^^^^^ J ^^^^^^ ^^^^^ ds + X, ^^^^^^ . 



"1 



(iMo) ■ 
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The deterministic homogenized solution uq is dominant in the Umit £ — > 0. As we will 
show, the remainder is 0{e) in L^{Q x [0, 1]). The term can be re-written 




1 1 1 



qe{s)G{x,s)ds, qe{s) :-- 



A{s) AoisY 

ns)-^){^,Io^^<it-l), 0<s<x 

Note that G{x, s) is piecewise continuous in s and so long as ""^ is bounded, G{x, s) 
is uniformly (in s) Lipschitz in x. We also show that Ve is 0{^/s) (in L^(Q x [0, 1])) and 
has a limit that can be characterized well. Note that these are slight generalizations of 
previous results. In particular this limit was studied by [8], and in the case of media 
with long-range correlations in [2] . This paper adds the feature that the random media 
is allowed to be non-stationary in one variable, and has no uniform (in u) ellipticity 
lower bound. 



3 Homogenization and Gaussian Corrector 

Here we obtain homogenization and Gaussian corrector results that are a generalization 
of known results (e.g. [17, 20, 7, 1]) to the case of media that has no uniform (in uj) 
upper or lower bounds. Proofs are deferred till section 7. 
First, we assume As{x) satisfies 

< z/i(w,e) < Ae(x) < Z/2(W,£) < DO. (8) 

We abuse notation by writing A^i^x) := A{x, |) = A{x, The form A{x, f ,w) 

emphasizes that A is a random field defined on a probability space (fi, S, P). The form 
A{x, ^,^,0) emphasizes the dependence on a special random vector ^ : ^ M"*, and a 
(possibly infinite) sequence of random variables 9. Once we fix ^, A(x, y) exhibits some 
stationary in y (although only weak-stationarity of A~^ is needed for homogenization 
in one dimension). 

Note that functions such as {A^} do not depend on y, so we write {Q} {x) = 
E^{Q{x,0)} for functions Q{x,^,y,6) whose conditional expectation does not depend 
on y. It is necessary to make some ergodicity assumptions on the process (in y). In one 
dimension, we require mean-ergodicity of E^ which we express through decay 

of the covariance (10). We also assume ||/||l2 < oo. 

3.1 Homogenization 

We assume weak stationarity of A"^ and ellipticity of ^o- In other words, we assume 

0<ci{0 <'^i{A-\x,y)} = E^{A-\x,0)} = Ao{x)-^ < C2 < oo, 
and that the conditional covariance 

Cov^(xi, X2, z) := E^ { [A-\xi,y) - E^ {A~^} (xi)] [A-\x2,y + z)-E^ {A'^} (xa)] } , 
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is independent of y G [0, oo). We also assume 

\CoY^{xi,X2,z)\ < r{z) with |r||ii =: C^-i < oo. (9) 
Note that this imphes 



/ / \Cov^{xi,X2,y-y)dydy\<eCA-i- 
Jo Jo 



This ahows 
Theorem 3.1. 



'E^\\Ue - Wo||l2([o,i]) < ^V^UWl^VCa^. 
Moreover, as e ^ 0, with ^ fixed, 

Us{x) — uq{x) ^ for every X, a.s. 
and if the ellipticity bound (8) is independent of s, 



(10) 



3.2 Gaussian corrector 

To quantify the rate at which the random coefficient dccorrclates, we introduce the 
following mixing condition (referred to in the literature as p mixing). 

Assumptions 3.1 (Mixing Condition). For two Borel sets A,BcM., let Sa, Sb be 
the sub sigma algebras generated by A{x,y), x & [0,1], y & A and y £ B respectively. 
We assume the existance of non-negative bounded and decreasing (/? : [0, oo) ^ R such 
that (p^/^ G that also satisfies 

sup{|E[r?„?76]| : % G G 5b, Er?^ = Erj^ = 1,E% = E% = 0} < (^(d(^,5)). 

After conditioning on ^ (e.g. fixing one realization of a(x,^,) if ^ = a + 6) we are 
able to partially characterize the limiting distribution of — uq. 

Theorem 3.2. If A{x, y) is stationary in y, and (si, S2) ^ Cov^{si, S2, 0) is continuous 
at {x,x), A~^ satisfies the mixing condition (3.1), and 

then 

Ueix) - uo{x) dist. f . 

— — ^-^ > v[x) := / G[x,t)a{t) dWt, 

Jo 

where Wt is a one- dimensional Brownian motion, G is given by (7), and 

/oo 
Cov^{t,t,q) dq. 
-00 
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Remark 3.1. The continuity on Cov^ is equivalent to mean-square continuity of x i->- 
A{x,0)~^ (see e.g. [18]). This means more-or-less that the media varies slowly with 
respect to x. This is where scale-separation comes in. 

Remark 3.2. The variance of the above expression is then given by 

C G{x,tfa'^{t)dt. (11) 
Jo 

Remark 3.3. The lack of a uniform (in u) ellipticity lower bound is dealt with through 
the bounds on sixth moments and the mixing conditions. 



4 Large Deviations 

Here we derive a rigorous "pointwise" large-deviations result for u^, and an approximate 
rate function. 



4.1 Moderate and large deviations for our problem 

For small enough e, the homogenized solution uq captures the bulk of the solution u^. 
The corrector attempts to capture some statistics of the term — uq = + Re- Our 
corrector result shows that 



Ue{x) - uq{x) 



^ C G{x,t)q{t,-) + 

Jo ^ 

^v{x):= [ G{x,t)a{t)dWt. 
Jo 



(12) 



By definition this means that, for any £ > 

Ue{x) - Uq{x) 



F^[v{x) > i]. 



A relevant question is whether or not 

P^[ue{x) > £] P^[uo{x) + y/ev{e) > £]. 



(13) 



(14) 



The inequality > £ (for i > uq) is a large deviation since (12) (or a law of large 
numbers result) shows that concentrates near uq. On the other hand, £'~^/^(u£ — 
Uq) > £ is a moderate deviation. Generally speaking, (14) does not hold. Instead, 
from (13) we can only rigorously infer something about moderate deviations, namely 

PdMx) > V^£] « PdM^) + V^vix) > 

For simplicity, we consider large deviations at only one fixed x G (0, 1) and often 
change the notation to uq, «£, (dropping the x dependence). 

Definition 4.1 (Rate functions). A rate junction / is a lower semicontinuous mapping 

(such that for all a G [0, oo), the sub-level set ^/(a) := {x : I(x) < q} is closed) 
/ : M" — [0, cx)] . A good rate function is a rate function for which all the sub-level sets 
^'/(a) are compact. 
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Definition 4.2. We say that a family of random vectors G satisfy a large 
deviations principle (LDP) with rate function / if for all F C M" 

- inf I{y) < liminf £logP[y£ G T] < limsupelogPfy^ G T] < - inf I{y). 

Above, r°, r denote the interior and closure of F. 

Considering now our solution G M, by choosing A = [i, oo) we obtain limiting 
upper and lower bounds on P^[us > £]■ We also often have, for i > limE^u^, 

lim e log P[ue >£] = -!{£). (15) 

We will find a rate function /^^ for Ue- We also find a rate function / for UQ + Ve, which 
proves to be more tractable. Simulations show that lu^ / for £ ~ 0{1) but not too 
large (see figures 5, 6, and 9). 

The rate-function gives the exponential rate of convergence in the sense that 
given £ G M, 5 > 0, there exists eq such that for e < eo 

P^K >^]< e-(^(^)-^)/^ (16) 

As an example, consider the corrector in our 1-d problem, ^/ev. It is Gaussian with 
mean zero and variance equal to eCc for some Cc (given by (11)). Therefore ~ 
+ v^f ~ J\f(uo,eCc)- Keeping the first term in an asymptotic expansion of the 
complementary error function, we have that for i > uq, 



V27r(£ — Mo) 



Thus, when 



we will have £logP^[Mo + v^'W > ^] ~ —{^ — uq)'^ /{2Cc)- Note that this is a "small e 
and large \£ — uo|" condition, as it should be. Comparing this to (16) we see that the 
Gaussian corrector captures the asymptotic tail behavior when 

(18) 

If (18) does not hold, the corrector cannot capture the tail behavior of u^. A cautionary 
note is in order here. The large deviations result captures the exponential rate of 
decay, and important algebraic factors in e are not captured. So, for finite e the rate 
function can be used for comparative purposes, not to estimate the true tail. It should 
be noted that so-called concentration inequalities provide a number of upper bounds 
on sums of random variables, often in pre-asymptotic regimes. Often these require 
fewer assumptions but do not claim to be tight. For example, Chemoff's bounding 
method is given in (20). See [6] for a survey, and [19] for an application to uncertainty 
quantification. For independent media, we derive rigorous upper bounds for finite e 
(e.g. (20)). For dependent media similar bounds arc available but arc more complex 
than the asymptotic bounds. Another reason for using asymptotic bounds is that a 
"fair" comparison such as (18) can be made. 
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4.2 Independent sums and basic definitions 

It is instructive to start here. Let X„ be random variables and define 

1 ^ 

n=l 

Now, for A > 0, 

P[Sn >i] = P[NSn > Ne] = e{i^^^>^,} < e-^^%{e^^^^'-} . (19) 
With vector-vahicd random variables in mind, we define 

Definition 4.3. The logarithmic moment generating function for a random variable 
Y is defined as 



A(F,A) = Ay(A): = logE{e^-^}. 



(19) allows us to conclude 

(20) 



^logP[SN>i]<-SUp 
A>0 



The above {Chemoff's) bound holds without any assumptions. Suppose the X„ are 
i.i.d. This leads to 



n=l 



and thus 



1 1 ^ 

-A{NSn,\N) = -5]logE{e^^i} = Ax(A). (21) 

n=l 

Note that Ax is convex since by Holder's inequality 

Ax(Ui + (1 - t)\2) = logE {(e^i^)*(e^^^)(i-*)} 

< log |e {e^i^}* E {A2X}(i-*) I = tAx(Ai) + (1 - t)Ax(A2). 

Definition 4.4. The Frenchel-Legendre transform of A{Y, •) is defined by 
A*{Y,i) = A^{i) : = sup [A • £ - Ay(A)] . 

Inserting (21) back into (20) (taking into account negative A) we have a large 
deviation upper bound with rate function A'^(-). The "trick" is to obtain a lower 
bound and thus show that this upper bound is tight in the limit N ^ oo. This indeed 
is the case and the result is 

Theorem 4.1 (Cramer, [13]). The sum Sn satisfies the LDP with good convex rate 
function A^^ (•). Moreover, (15) holds. 

Since Ax(0) = 0, we always have A'^{£) > 0, i.e. we never have exponential 

growth. Jensen's inequality shows that Ax (A) > XEX, so for i = E,X we also have 
X£ — Ax (A) < 0. Therefore A^^ (EX) = 0. This makes sense in view of the law of large 
numbers. 
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4.3 The Gartner-Ellis theorem and contraction principle 

Definition 4.5. For any function H taking values in (—00,00], we define T>h '■= {x : 
H{x) < 00}. 

The Gartner-Ellis theorem will be used to prove results for our model problems. 

Definition 4.6. Convex A : ^ (—00,00] is called essentially smooth if T>\ is 
nonempty, A is differentiable throughout P^, and A is steep. Steep means that 
limj_^oo |VA(Aj)| = 00 whenever {A^} is a sequence in converging to a boundary 
point of V\. 

Theorem 4.2 (Gartner-Ellis [13]). Suppose 

A(A) : = limelogEje^"'^-^^} 

exists as an extended real number. Furthermore suppose that A is essentially smooth, 

lower semicontinuous and that the origin belongs to the interior of T>\. Then 
satisfies an LDP with good convex rate function A* defined by 

A*{£) : = sup [A-^- A(A)]. 



Remark. If n = 1 and A* is finite in a neighborhood of ^ > a := limEZ^ then convex 
A* is non-decreasing on [a, 00) and continuous in this neighborhood. Therefore (15) 
holds. 

Notice that the Gartner-Ellis theorem does not require independence. Indeed, one 
can use it to prove an LDP for mixing random variables [13, 9]. 

While the Gartner-Ellis theorem allows us to obtain an LDP for oscillatory integrals, 
we need the contraction principle for functions of those integrals such as 5. 

Theorem 4.3 (Contraction principle). Suppose f : M" — t- M™" is continuous and I : 
[0, 00] is a good rate function for the family of random variables Z^ and associated 
measures /x^ (HeiA) = P\Z^ G J^). For y G W^, define 

I'{y) ■.=\ni{I{x) : x^W\y = f{x)]. 

Then I' is a good rate function controlling the LDP associated with the measures //go/"-'- 
(^ieof-\B) = P[f{Z,)^B]). 

In other words, with := f{Z^), the rate at which Z^ concentrates away from £ 
will be determined by the point in f''^{{i}) holding the most mass. 

4.4 Large deviations for 

In light of (5), we have 

Ue{x)=g{Ze), : X (0, 00)^ M, g{z) = g{zi,Z2,Z2„ZA) = -zi + Z2 — , 

z^ 

Z, = (^,i,...,Z,4), Z,, = I'^ds, 

Jo Ae{s) 

i?l = Fl(o,^), H2 = F, if3 = l(0,x), ^^4 = 1. 

(22) 
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Also set 



H : — {Hi,H2, H^, H4), 

so that X- Ze = fx- H/Ae. 

Using the Gartner-EUis theorem we wih find rate function Iz^ for Z^, and then by 
the contraction principle the rate function for is 

IuM= ^^i,JzM- (23) 

The real work involved here is in characterizing the limiting Cramer functional 
lime log exp{e~^A • Z^}. We will do this for the media types (i) and (ii) mentioned 
in the introduction. 

4.4.1 Pcirameterized, independent, uniformly elliptic media 

We assume here that the high frequency media is piecewise constant and independent: 

A{x,y) = a{x,$,) + b{y,9), with every realization a{x,^) continuous, 
00 

Hy,G) = nJ2^^Mn-iMy)^ ^^-i.i.d. 1^1 (24) 

0<iyi< A{x,y) < V2. 

Considered as a discrete process (at points centered at ne), the field is stationary 
and ergodic. It is not truly stationary since the correlation W.^A{x, y)A{x, z) depends 
on more than the difference y — z. Nonetheless we apply our theorems and obtain 
results that are validated by simulation. The correlation condition is satisfied trivially 
since A{xi,yi), and A{x2,y2) are conditionally independent whenever \yi — y2\ > 1. 
Note also that the low and high frequency parts are in the form of Karhunen-Loeve 
expansions, but the total field is not. 

Since this solution (5) involves A^"^, it is not surprising that we need to define 

■■ = -^-E ~ T^vM = ( - ( - - a ) ) , (25) 

The main result for this media is 
Theorem 4.4. With g, Z^, H(s), given by (22), (24), (25), define 

A(A):= [\iVa^s,O^X-Ji{s))ds, 
Jo 

then A G C°°(M^) is a convex function such that when A^ is defined by (24) 

lime log E^e""'^-^^ = A(A). 

Moreover, for fixed Z^ satisfies a large deviation principle with good convex rate 
function 

A*{£) : = sup [A-^- A(A)], 
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and satisfies a large deviation principle with good rate function 



Proof. We will condition on and then approximate a{x) by a piecewise constant 
function so that a{x) appears only as a parameter. Now 

So for AT < £-1 < AT + 1 G N 

1 ^ 

e log |e^"'^-^^ } = eA^— ^ log |e^-^^(") } + e log E^ je^'^'^ } , 

n=l 

H(A) z-^- H(A) 



Since i/^^ < (a+fb^)"^ < i^f ^ we also have upper and lower bounds on E^ {exp {A • Yn}}- 
Therefore elogE^ {exp {A • Y^}} 0. Using also the fact eN ^ 1 we have 

1 ^ 

lim£logE.|e^"'^-^4 = lim — VlogE. |e^-^^(")]> . 



We henceforth study the limit on the right. 
We can approximate 

g Q 

Hiv(n) : = max H(— ), ajv(n) := min a(— ). 

n— l<s<n iV n— l<s<n iV 

We thus have 

1 ^ _ 
£logE{Ai;j£} < - ^A(K^(„),A•H^^(n)). (26) 

n=l 

Similarly we can choose HjY(n) and aAr(n) to provide a lower bound. Together they 
yield 

N , W 



n=l n=l 

1 ^ — 



n=l 

The ellipticity bounds imply that (a, A) A(V^,A) is C°° for a taking values in 
the closure of Range(a). Therefore, 



A(F„^(„), A • HAr(n)) ^ A(K( « ), A • H(-)), ^ oo. 
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And thus (after extending a^, Hjv to be piecewise constant), and using the continuity 
of s !->■ a(s,^) and H(s) we have 

1 ^ _ /.I _ 

JaT.^('^-Nin)^^-tiNin))= / A(y„^(,^),A-H^(siV))ds 

n=l 

^ rA(K(,),A-H(s))ds. 
Jo 

The same holds for the lower bound. We have thus shown 

hm e log e'~'^-^^ = C A(F„(,) , A • H(s)) ds. (27) 

JO 

Since E {exp {AV^}} is finite for all a, the hypothesis of the Gartner-Ellis theorem 
are trivially satisfied. Recalling the definition of V^, A, we see that the theorem is 
proved. □ 



4.4.2 Convolved media with no uniform (in uj) lower bound 

Here wc obtain a large deviation principle for dependent media given by a convolution of 
random variables that, while being positive, have no uniform lower bound. Convolution 
provides a convenient way to generate dependencies. 

To avoid additional technicalities, we restrict our attention to families of functions 
indexed by e such that £ N. Starting from this result, using the notion of exponen- 
tial equivalence, and adding an assumption of a finite logarithmic moment generating 
function, it would be possible to prove an LDP for general e £ (0,oo). 

Define 

^ oo 

A^{x) = A{-), where -j-- := ^ l[n-i,n)(s)7n, 

1 ' (28) 

ln-= X] hn-mPm, K > 0, := ^ /ife < OO, 

rra=— oo k 

and the {/3m}m=-oo non-negative i.i.d. random variables each depending on the 
same random vector (parameter) ^. The random variables 7„ are well defined so long 
as the characteristic function (f)M{t) ■= n|fc|<Af ^ jg^^'^fe/^j has a continuous hmit (/)(t). 
We take this as an assumption and proceed. 

In this case, the large deviation principles for Z^, direct result of lemmas 

4.1, 4.2 (below), the Gartner-Ellis theorem, and the contraction principle. 

Theorem 4.5. With g, Z^, and H(s) given by (22), and (24), define 

A(A):= rA^(^)(||/i||iA.H(.))rf., 
Jo 

then A G C°°{M.^) is a convex function such that when A^ is defined by (28) 

lime log E^e''"'^'^^ = A(A). 
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// in addition A is steep (see proposition 4-V' then for fixed ^, satisfies a large 
deviation principle with good convex rate function 

A*{i) := sup [A-^-A(A)], 
and satisfies a large deviation principle with good rate function 



Lemma 4.1. Let 'jn := Y^j^Ki-kl^k with non-negative {/ifc} G £^ and non-negative Pk 
i.i.d.with the same law as p. For G G L°°([0, 1]), define 



Sn •■= G(s)7uArj 

Jo 

Then for A G M the following limit exists in (—00, 00] 
limlA(iV5jv,A) =^'a^(AG( 

with A/3 as in definition 4-3. 



ds. 



1) ds, 



Proof of lemma 4-1- For a result involving sums of convolved stationary random vari- 
ables see [11]. The key difference is that here we allow the moment generating function 
to be infinite, and the term G{s) makes the sum non-stationary. This necessitates a 
new proof technique that relies on non-negativity of the P^ a^^d the h}~. 



NSn = N G(s)7L3ivj ds = V / G{s/N) di 

Jo Un-l 

= Y, / G{s/N)ds ^Pkhn-k. 



Making the substitution j = n — k we have 



N-k 



NSn = YPkHN,k, HN,k := X? ['^^ G{s/N)ds. 
fcez j=i-k Ji+k-i 

This is a sum of independent random variables. Since P^ ~ P, 

11 1 /"°° 

-A{NSnA) = j-J2M^H^,k) = 1^ / A0{XH^^^,i)ds. 

Using (43) and changing s Ns we have 

—A{NSn,X) = J (^jv(s)ds, (^jv(s) := A^ (AiJjv, [sJVj ) ■ 
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The key to our proof is the fact that i?Ar,fe itself can be written as an expectation 
HN,k = / Gis/N)h^,_k+ii = E,, {||/i||iG(-/iV)l[o,iv]} , 

^0 



fis) 



h 



[s-k+l\ 



ds. 



Since vr^ is a density that concentrates near s = k — 1, for almost every s G [0, 1], 
'^T^lsN] {ll^lli^(V^)l[o,Af]} ~^ ll^lliG'(s) (this is a result on the Lebesgue set for inte- 
grable functions, see e.g. [15] theorem 3.20 or [13] theorem C.13). Also, 



\Hn,IsN\ \ < ll<-^l|oo ^ hj. 

j=ll-sNi 



(29) 



(30) 



Since this tends to zero for s ^ [0, 1], we have 

To prove the lemma, we thus have to show / (fN J ^■ 

We first obtain an upper bound in cases where A(A) < oo. When A(A) = oo, the 
lower bound we derive will be infinite, and thus the upper bound as well. 

^N{s)ds = / A^ (e^^^^j {A||/i||iG(-/iV)l[o,iv]}) ds 

-oo J —oo 



-oo 
f oo 



/oo 
E^^^.j {A^ (A||/i||iG(./iV))} ds 
-oo 



TV 



(A||/i||iG(t/iV))7rL,Arj(t)di 



-1 /-co /"iV 

-/ / Afi{X\\h\UGm))^[si{t)dt 

%j — oo 

/oo r f'\ 
/ A^(A||/i||iG(i))7rL,j(tAr)di ds 
-oo UO 



ds 
ds 



(31) 



= / Kp{\\\h\\,G{t)) 



■Jo 



[ Trisi{tN)ds) dt 

J —oo J 



f\^{X\\h\\iG{t))dt= [\{t)dt. 
Jo Jo 



Above, the inequality is due to Jensen's inequality and the convexity of A^. The change 
of integration order is justified by Fubini since ip ^ L?-. 

We now obtain a lower bound. We will show that (pN > ip, with ip ^ L?-. Then a 
corollary of Fatou's lemma gives us 

liminf j (/?jv(s)ds> j (^(s) ds = A(A). (32) 

With M(t) := E {e*^}, we note that 

^n{s) = \ogM{\HNXsN\) > logM(-|Ai7^,L-jvj|) = - log l/M(-|Ai7jv,L.ivj I)- (33) 
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Also, using M(— |Aii"jv,[sArj I) < 1 ^i^d (29) we have 

log i/M(-iAi^.,,.j I) ^ log (i + 



^ |l-M(-|Aff;v,L.ivj|)| ^ |1 - M(-|AF;v,L.ivj 



E|i-e-'^^^>L»ivjl| ^ |Ai/jv,Uivjl ^^^^ 



< ' . , , , ' < 



So, 



- e^;3(-A||/i||i||G||oo) - e^/3(-^ll'^llillG'l|oo) 

[N+l-sN] 

<e-AM-^IWi||G|k)||Gr||^ ^ h^, 

j=ll-sN\ 



logl/Af(-iAi7^,L«^j|) < e-^''(-^ll'^llill^ll-)||G||oo||/t||i. (35) 



This bound works for all s. However, when s is away from [0, 1] the summation is over 
the tails of hj and we can do better. Specifically there exists A^O; -Ai such that for 
N > Ni> No and s [-1,2] 

[N+l-sN] [No+l-sNoi 
j=\l-sN\ j=Ll-sAroJ 

Therefore, using (34), for s ^ [-1, 2], AT > ATi, 

\_No+l-sNo\ 

logl/M(-|AiJjv,L.ivj|) < e-^''(-^ll'^llill«ll-)||G||oo (^6) 

3=\l-sNo\ 

So define 

.e-M-M\hh\\G\\^)\\G\\^\\h\\i, se [-1,2] 

^INo+l-sNo 
^j=ll-sNo\ 



^^'^ - ~ -e-^^(-^IM^II«ll-)||G|UESri;;^^^ h„ s i [-1,2]. 



Then for A > A^i, (33), (35), and (36) show that <y?Ar > V- To show G J-^ , we note 
that 

. nOO L^o + i-siVoJ LAf() + l-sJ 

/ ms)\ds<C J2 ^id. = — / Y hj'^' 

Js^l-1,2] J-^ ,=lI^7VoJ "^-^ i=Li-.J 

= C\\h\\i <oo. 

We thus obtain (32) and the proof is complete. □ 
Lemma 4.1 allows us to obtain a limiting Cramer functional for the convolved media. 
Lemma 4.2. With H defined by (22) set 

A(A):= rA;3(||^||iA-H(s))ds 
Jo 

Then 
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(i) Restricting attention to e such that e ^ eN, 

lime log Efe^"'^-^^ = A(A). 

(ii) A is lower semicontinuous. 

(Hi) Let b G (—00,00] be the number satisfying {A : Eexp{A;S} < 00} = (—00,6) or 
= (—00,6], and let Hx be the value o/H(s) that maximizes A • H(s). Then 

{A : ll^lliA • H;, < 6} C Pa C {A : \\h\\iX ■ < b}, 

with A' finite throughout . 
(iv) If A is steep, then A is essentially-smooth. 

Proof of lemma Write 

elogE^e"^"'^-^^ : = elogE^ |exp \e~^ "^^^^ dsjj . 

Applying lemma 4.1 with A = 1 and G = A • H we obtain that A is indeed the form for 
the limiting Cramer functional, thus proving (i). 

To prove lower semi-continuity, we note that whenever A„ — >■ A and A(A„) < a then 

a > lim inf A(A„) = lim inf / A^(||^||iA„ • H(s)) ds 

Jo 

> [ A^(||/i||iA-H(s)) = A(A), 
Jo 

by a corollary of Fatou's lemma (using the fact that the integrand is bounded below). 

To show (iii), suppose first ||/i||iA • < b. Then A^(||/i||iA • H(s)) is bounded and 
differentiablc and therefore A(A) is too. If on the other hand b < \\h\\iX ■ Hx, then 
since H is continuous, A^(||/i||iA • H(s)) equals +00 on a set of positive measure, hence 
A(A) = 00. We thus have our bounds on V\ and it follows that P^ is non-empty and 
A' exists in P^. 

Lastly, (iv) follows from (ii), the assumption of steepness, (iii), and the definition 
of essential smoothness. 

□ 

As hinted at by lemma 4.2, steepness of A is a condition that needs extra work 
to check. We formulate a necessary and sufficient condition below, and then three 
sufficient conditions that are easy to check. 

Proposition 4.1 (Steepness criteria). Extend A'^ to map R [0, 00] by setting A^(i) = 
00 whenever A^ {t) = 00. Then define ifj : — >■ R 6?/ 

/•X i-l 

Ki{r]) : = / A^(F(s)r?i + 7/2) ds, K^iv) := / A^(i^(s)?7i + m) ds. 

Jo Jx 

Then A defined in lemma 4-2 is steep if and only if Ap is steep and for every rj G dV^,, 
Ki{rj) = 00, i = 1,2. 

Moreover, with Fm '■= max^ F{.s), F^ ■= min^ F(s), A is steep whenever one of the 
following sufficient conditions hold: 
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1. Ap is finite everywhere 

2. F is piecewise and = (—00,6). 

3. Let {si, . . . , Sn} be the points where F{si) = Fm or F{si) = Fm- Then there exist 
neighborhoods Ni = (sj — 5, Sj + (5) fl {s : F^ < F{s) < Fm} such that on Ni, F 
admits an expansion of the form 

F{s) : = F{si) + c{s - SiY + R{s - Si), 
R(s-Si) ^ R'(s-s,) 

7^ — ^ ^ 0, , ^ . ^0, N,3s^ Si. 

[s-SiY [s-SiY 1 

Then with r := min{ri, . . . , r„}, A, b = dV^^, A is steep if 

Kit) 



L 



Remark 4.1. Since ^ > 0, we have A'^ > 0, and therefore Ki : [0, 00] are well defined 
as integrals of functions taking values in [0, 00]. 

Proof. We note that the first condition "A^ is finite everywhere" trivially implies that 
A is steep. If A^ is not steep, then it is easy to construct an example showing that A 
is not steep either. So from now on we assume A^ is steep but 7^ ^■ 

We now show the necessary and sufficient condition involving the Ki. Recall 

H : = (Fl(o,j;), F, l(o,x), !)■ 

To show steepness we must fix A G let 9 A" ^ A and show |VA(A")| 00. 
Define the function T : ^ [0, 00] by 

r(A):= /\^(||/.||iA-H(.))d.. 
Jo 

We claim that A is steep if and only if r(A) = 00 for all A G dV^. Indeed, if r(A) < 00 
for some A G dV^, then since for t G (0, 1) t i->- A^(||/t||iiA • H(s)) is finite and non- 
decreasing (by convexity), 

|VA(a)| < |||H|||Lcor(a) < |||H|||Loor(A) < oo, 

and therefore A is not steep. On the other hand suppose r(A) = 00 for all boundary 
points A, then choose one along with a sequence 5 A" — ^ A, then by Fatou 

liminf |VA(A")| > liminfr(A") > r(A) = 00. 

Therefore A is steep. 

We now show that T>^ = V\. Note that A G implies A^(||/i||iA-H(s)) is bounded 
for all s G (0, 1), and AG (A — A + 5) for some (5 > 0. Since = T)°^, , the same 

holds for Ap{\\h\\i\ ■ H(s)). Hence C V\. A similar argument shows V\ C V^. 

We now have A is steep if and only if r(A) = 00 for all A G dT)\ = &D^. Our 
next step is to change variables to simplify this boundary. To that end, note that 
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r(A) = K,{\\hUX, + A2), ||/i||i(A3 + A4)) + K2i\\hUX2), ||/i||i(A4)), so by a change of 
variables r) = (Ai + A2, A3 + A4, A2, A4) we have 

r(A) = i^i(||/i||ir/i, \\h\\iri2) + i^2(||/i||ir/3, WHiVa)- 

Changing variables again ry i-)- ?7/||/i||i and taking note of the non-negativitiy of the Ki, 
we see that A is steep if and only if Ki{rj) = 00 for all 9 77 € dV^i- 

Having proved the necessary and sufficient condition, we use this to show the three 
sufficient conditions. The first has already been shown. 

As for the third sufficient condition, choose an extremal point s G (0, x) and assume 
rj € T^Ki ■ Without loss of generality, assume we have the given expansion in the open 
set {s — 6, s). Now Ki{r]) = 00 if for all 6 > 



A'g{F{s)rii + m)ds = <x. (37) 



We shall reduce this condition to the type stated in the proposition. 
Our expansion gives us 

t{s) : = F{s)r]i + r]2 = b - c\s - sr + ^(s - s), 

for a new positive constant c and a new function R (differing from the old i? by a 
constant). After possibly shrinking 5 we can solve for s{t). 

s{t) = s- c-^''{h + R{s - s) - tf'\ 



Therefore, (37) holds if and only if 

r A^(F(s)r7i+r72)ds= f Mg{t)s' {t)dt = rx^. (38) 

Js-5 Jt{s-5) 

Differentiating we have 

^ ' V cV'-(6 + R^s - s) - t)(^-i)/^ ) c^^ib - t)('-i)A (6 + R{s - s) - t)(^-i)/r ' 

Noting that b + R{s — s) — t = c{s — s)'^', and using our hypothesis on R, we have 
positive Cl , C2 such that 



(6 - t)('^-i)A - ' ^ - (6 - t)(''-i)A ■ 
Due to (39), (38) is equivalent to 



dt = 00. 



(6-i)('-i)A 

Sufficient condition three then follows by considering all possible such points s. 

Sufficient condition 2 follows from 3 since given 2 we have the expansion in 3 with 
r > 1 and therefore 

/ 77 w iw dt< / A'.(i)di = lim / AUt)dt 

= limA^(,5)-A;3(6-l) = oo, 

in light of our assumption = (~oo,6) and the lower-semicontinuity of A^, which 
follows from Fatou's lemma. □ 
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5 LDP approximations and generalizations 



5.1 Rate functions for approximate solutions 

To compute (naively) the rate functions for requires first a four-dimensional (convex) 
optimization to obtain A*, and then another four dimensional optimization to obtain 
lue ■ Our corrector theory shows that one can rigorously approximate = uq + in 
an ^/£ neighborhood of uq. Motivated by this we consider the rate function for uq + v^- 
However, a large deviation will necessarily take us outside the ^/e neighborhood, so 
more discussion is in order. Looking at terms in the expansion (6) = uq + + Re 
we see that when terms of the form 

1 r^i^i'i,. 



(Hk/Ao) Jo Ae{s 

are not too large, we can approximate ^ uq + Vs- An exact rate function for uq + Vs 
can then be calculated. We call this our approximate rate function Note that this 
is equivalent to approximating g (from (22)) by some map g and using the contraction 
principle. It could be argued then that the inverse images g~^{f-} ~ g^^{i} when 
i is sufficiently small (or some other better conditions). Then continuity of the rate 
function A*{z) shows that lu^ii) « Since we can also represent Ve{x) (for fixed x) 
in terms of an integral against a single function G{x, s), we can obtain a rate function 
for no + Vs without using the contraction principle. This gives us a more explicit form, 
and since rate functions are unique (lemma 4.4.1 [13]) these methods give the same 
result. 

We present now the LDP for uq + Vs- The proof is a simplified version of the LDP 
proof for Ue- 

Proposition 5.1. With G, given by (7), we have: 

lime log E^e^"'^("o+''^) = A(A). 
Where, when is the parameterized media defined by (24), and Va given by (25), 

A(A): = Ano+ r [-A^ + A(V;(,),AG(s))l ds, 
Jo L ^o(sj 

and when A^ is the convolved media defined by (28), 

A(A) : = Ano + [' f-A-^ ds + A^(||/i||iAG(s)) 



ds. 



In either case A is convex and whenever A is steep uo + satisfies a large deviation 
principle with good convex rate function 

: = sup [A^ - A(A)] . 

Moreover, whenever I is finite in a neighborhood of i > limEu^ then (15) holds. 

Remark 5.1. The steepness criteria in proposition 4.1 also apply here with G{x,s) 
(written G{s) when we fix x) replacing F{s). 



20 



5.2 Generalizations to incompletely characterized media 



Gaussian corrector results, e.g. theorem 3.2, require knowledge of the first two moments 
of the random media (along with other "niceness" assumptions such as mixing and 
bounds on higher moments). The question arises: How much must be known about 
the random media for a large deviation result? Here we partially answer this question 
and leave a thread open for future work. 

LDP results for mixing random variables are available [13, 9]. These give existence 
but not an explicit form for limiting Cramer functional A, and hence only existence 
of an LDP (instead of an explicit form). Restricting our attention to the case of 
parameterized media (section 4.4.1) or convolved media (section 4.4.2), we ask, "how 
general can the assumptions on the 9^ or (3^ be?" In general, we cannot expect the 
moment generating function of our media to be known. Notice that 



A > A ^ sup 



M-k{\) < sup [A£ - A(A)] . (40) 



In other words, if we can obtain upper bounds for the moment generating function, 
then we can obtain a lower bound on the rate function. So a strategy would be: First, 
write the limiting logarithmic moment generating function A(A) := lim£_>.o eA(e~^i'£, A) 
in terms of the logarithmic moment generating function of the random media (e.g. A^ 
in lemma 4.2). Second, find upper bounds for A^ using e.g. Bennett's inequality 
(Lemma 2.4.1 in [13]) (this bounds the moment generating function of a bounded 
random variable in terms of its mean and variance). Third, a rigorous limiting upper 
bound is now available via (40). 



6 Model Problems 

Here we explore two specific examples and give numerical results. In both cases the 
right hand side f{x) = 1 for a; G (0.45,0.55) and f{x) = elsewhere. Hence F{s) 
is piecewise smooth. Therefore, using sufficient steepness condition 3 (proposition 
4.1) our logarithmic moment generating functions will be steep (since s i-> G{x,s) is 
piecewise smooth). 



6.1 Numerical results for parameterized media 

Here wc use a field that fits into the framework of section 4.4.1. This gives some control 
over the large deviations. 

Let ^ = (^0) ■ ■ ■ 5 ^7) be ~j J W[— 1, 1]. We then set 

a{x, = max 1 1 + 2 ^ ~ ^ Ur"^ sin[(2m + l)Trx], ^| . 

l ' m=0 ) 

Next, let 6 = {9i, . . . ,0m, ■ ■ ■) be an infinite collection of identically distributed 
independent rvs 9i ^ U[—l, 1] (which are also independent of the ^i). Put 



^ 00 



2 

m=0 
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In other words, b{y, 6) = 9^ when m < y < m + 1. 
We are ensured of the ehipticity condition 

1 . , X. 7 
— < A(x, -) < -. 
32 - ^ ' e ^ - 2 

The resultant media is pictured in figure 1. 




Figure 1: Comparison of two realizations of parameterized diffusion coefficients (section 6.1). 
The homogenized coefficient differs significantly only when a{x, ^) is small. In both cases, 
e = 1/50. On the left a{x,C,) is such that no values of 6i ~ 1, 1] bring close to zero. 
This realization of a(x, ^) is quite typical and is referred to as our "mild" coefficient. On the 
right we have a less common "wild" coefficient. 



Following as in section 4.4.1 we characterize Va := (a + 0/2)-^ with ~ U[-l, 1]. 
We have an explicit density for Va, 



1 



^1/2 



<Va<- 



-1/2 



Hence, 



Therefore, 



E 



1/4- 



I ' r' 




\a + b{0) J 





log 



a+ 1/2 
a - 1/2 



Cova(r) 



0, otherwise. 



(41) 



These calculations are enough to give us explicit integrals defining the homogenized 
term uq and the corrector v,,. Namely, uq solves 

^Ao{x)^uo = /(x), 
ax ax 
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and the corrector is given by theorem 3.2 with a'^{t) as above. 

The large deviations result for is given by theorem 4.4, and our approximate rate 
function by 5.1. In particular, the limiting Cramer functional (theorem 4.4) is given 

by 

„1 .(a(s)-l/2)-l ^X-Uis)v 

A(A) : = / log / 2 — d-uds, 

Jo J(a(s)+l/2)-i V 

and for the approximate rate-function, 

A(A) : = Auo + / -A-r-H + log / 9— ds. 

Jo [ Aq{s) 7(a(s)+l/2)-i 

Wc define empirical rate functions £{u£{^.b),t), £{uQ{Q.b) + ^/ev{Q.b),^) (for ^^(O.S) 
and iio(0.5) + y^f (0.5) respectively) as follows. With {Xi, . . . ,X7v} a set of samples 
from W {W = Us OY W = uq + s/sv), we set 



_ 1 ^ 

W:= lim — VXj-, 



. _ J ^ log F Ef=i lx>K^i), ^ > w 



S{W,i) : 



(42) 



Note that limTv-^-oo ^(■f^o + y/£v,£) has an explicit expression, and we use this in place 
of (42) to compute £{uq + y/£v,i). Since Xj > £ is a rare event wc cannot compute 
£ by direct sampling. A crude importance sampling technique was used whereby the 
6i ^U[—l,l] were replaced by with (scaled and shifted) Bradford random variates, 

""^^-"^^^ = 21og(l + c)(l + f(^? + l))'l-l<^^^^- 

As c > increases, the draws 9i are more likely to concentrate near —1. This gives a 
smaller diffusion coefficient and hence larger solution. Calculation of £ must then be 
re-weighted by the factor Truniformi&i) /T^Bradi^i)- See e.g. [12, 21] for an overview of 
importance sampling. 

In figure 1 two realizations of the parameterized media are shown. We will fix 
the low frequency part a(x,^) and study the behavior of the solution over different 
realizations of 6(|,0). The "mild" medium (left) has a(x,^) far from zero, so no 
matter what b{j,9) is the solution is small. The "wild" medium has a section of very 
small a(x, ^). In all cases, one notes that the homogenized coefficient Aq(x,$,) differs 
from the low-frequency coefficient a{x, ^) most when the medium is small. In this case, 
Aq < a in an attempt to affect a large jump in the solution uq to approximate the 
often large solution (although b{y, 9) is symmetric about 0, the resultant solution 
Us is not symmetric about uo). 

We verify theorem 3.2 in figures 3 and 4. In figure 3 one can see that the pdf of 
the corrected solution (at the fixed point x = 1/2) agrees well with the true solution 
so long as e is small enough. The fit is worse for the "wild" medium, and in particular 
the true pdf shows an asymmetry that the Gaussian corrector cannot have. In figure 
4, one sees that when e ^ 1/50 (or smaller) the corrector captures the variance quite 
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l/£ X 



Figure 2: Left: convergence of E{||Me — uo||} and E — where Ua is the result 
of a truncated coefficient expansion (using a{x,S,) only). This verifies theorem 3.1. Right: 
Homogenized solution mo(x, ^) and many realization of the corrected solution mo(x,^) + 
y/ev{x,uj). For all realizations ^ was fixed at the same value as the wild coefficient (figure 1 
right) 




Figure 3: Left: Observed density of ^(0.5) for the mild medium (figure 1 left). Plot shows 
that the corrector captures the bulk of the variance quite well. Right: The wild medium 
(figure 1 right) is shown. Results are not as good. 
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Figure 4: Verification of tlieorem 3.2. Conditional variance {^^(a;, w)^} and 
E^{(uo(x) + ^/ev{x,UJ))'^} for various e. For all realizations ^ was fixed at the same value 
as the wild coefficient (figure 1 left). Plots show good agreement when 1/e > 50. Note also 
that {{uo{x) + y/ev{x)y} is available explicitly via (11). 

well. Rate functions for the mild medium are compared in figure 5. One can see 
that the corrector and true rate function are almost indistinguishable until the true 
solution saturates around 0.35. The approximate rate function also works well up until 
ti(0.5) ~ 0.3. In this case one could use the approximate rate function to see a-priori 
that the corrector stands a chance of capturing the large-deviation behavior well. The 
case is different for the wild medium LDP results in figure 6. Here one can see that 
the corrector rate function separates from the true rate function fairly early on. While 
the fit between the approximate rate function I and the true rate function £{u£, •) is 
not perfect, one could still tell, using only /, that the Gaussian corrector stands little 
chance of capturing the large deviation behavior. 

It should be noted that since for the mild medium, the maximum possible value 
of ^^(O.S) was approximately 0.35, we consider i = 0.03 a large deviation. The scale 
is harder to set with the wild medium since our sampling could not achieve results 
near the maximum. However, one does note (figure 6) that by £ ~ 0.6 the empirical 
rate function differs from a Gaussian rate function by quite a bit. For that reason, 
we consider i > 0.6 to be a large deviation. It is important to note however that the 
approximate rate function, being based on a linearization, does differ from the true 
rate function for large enough i. 
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Figure 5: Left: Comparison of empirical rate function for m^, S{us,i), with the corrector 
rate function {i — uo)^/{2Cc) in the case of the mild medium (figure 1 left). In this case, 
the corrector captures the large deviation behavior well up until ~ 0.3 when the true solu- 
tion nears its theoretical upper bound. Right: Comparison of theoretical approximate rate 
function with empirical rate function S with the mild medium and e = 1/100. The 
approximate rate function works quite well for values ~ 0.30. Since the largest possible 
solution is ^^(O.S) ~ 0.35 we consider 0.03 ~ 0(1). 
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Figure 6: Left: Comparison of empirical rate function for u^, £{uei i), with the corrector rate 
function S{uq + y/6v,i) in the case of the wild medium (figure 1 right). Right: Comparison 
of theoretical approximate rate function with empirical rate function S{us,i) with the 
wild medium and e = 1/100. 
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6.2 Numerical results for convolved media 



Here we implement a particular case of the media described in section 4.4.2. With 
e^^ G N, we define 

^ oo 

A{x) = A{-), where := ^ l[„-i,n)(s)7n, 



n=l 



In ■ ^ ^ hn—mf^mi 

m=—oo 

hn>0, ||/i||i := ^/ifc < OO, 



and the {/3m}m=-oo i-i-^- chi-squared random variables with ^ degrees of freedom. 
This means 



/3i~vr^(/3|0 



^5/2-lg-/3/2 



r(^/2)2«/2 ■ 

The moment generating and characteristic functions of every /3„ are 




0.2 0.4 



0.8 1.0 




Figure 7: Typical realization of random media when = 1 (left) and when = 3 (right). 



H / (l-2A)«/2' H / (1-2; 



(1 -2it)€/2' 



(43) 



The random variables 7„ are well defined since the characteristic function (j)M{t) 
Y[\k\<M'^ {^^^^^^} ^ continuous limit (j){t). Indeed, 



^ \og{l-2ithk). 

k\<M 



log (j)M{t) = - 

This converges absolutely as can be seen using | log(l— 2zt/iyfc)| < C\2thk\ and \\h\\i < oo. 
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Note that E^/S^ = ^, E^{pn - 0^ = 2^, so 



= 2gy^/t_fc/i|^|_fc. 



So the single random variable ^ defines the coarse-scale randomness. From realization 
to realization ^ varies with a geometric distribution (with parameter 1/5) i.e. 

1 /4^f-' 



Truncation has no meaning in this context, so we consider homogenization. We 
have 

A^^ = E {Ae^'^} (x) = ^\\h\\i = const., (45) 
and then uq is the solution to 

d^ 

-C\\h\\i^uo = fix), no(0) = txo(l) = 0. 

Using (44) we obtain o"^(t) = 2^||/i||^. The Gaussian corrector is then given by 
theorem 3.2. The large deviations result and rate function is given by theorem 4.5 




Figure 8: Pdf of vs. that of uq + ^/6v for two values of e and the convolved media. On 
the left, K = 1, and on the right k, = 10. In all cases ^ = 1. Shows good agreement once e is 
small enough, although agreement is worse when k = 10. 



with 



A(A): 



2^°sr 



lA-H(s) 



ds 
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and for the approximate LDP (proposition 5.1) 



A(A) : = Alio + 



ds. 



Figure 7 shows typical realizations of the diffusion coefficient when the media "build- 
ing block" (3 has ^ = 1 or ^ = 3 degrees of freedom. More degrees of freedom means 
larger Us- The behavior however is not analogous to the "mild/wild" comparison of 
section 6.1. In particular, the corrector captures the bulk of the distribution (moderate 
deviations) for all values of ^ so long as e is small enough. For this reason we only 
picture pdfs for = 1 (figure 8). This is expected since a x| random variable behaves 
similar to a Gaussian random variable M{(,, 2^) when is large. The media correlation 
length K does effect the results, and in particular k must be much smaller than 
for the corrector theory to work well. This is expected since sums of highly correlated 
random variables tend to a Gaussian at a slower rate than independent ones. Figures 9 




Figure 9: LDP results for the convolved media. Comparison of empirical rate function 
S{ui;,i) with theoretical (approximate) rate function Here ^ = 1, \\h\\i = 1, n = 1 and 
e = 1/100 (left) or e = 1/10 (right). 

and 10 show that the approximate rate function captures the large-deviation behavior 
well (for i not too large). As in the case of the wild medium from section 6.1 we con- 
sider i a large deviation if the empirical rate function £{u£, •) differs significantly from 
the Gaussian rate function at that point. In all cases, e must be small enough, but 
when it is the match is good. Figure 11 shows that the corrector cannot capture the 
large deviation behavior of this media. As is the case with the parameterized media 
of section 6.1 the approximate rate function does a much better job than the corrector 
(once 8 is small enough). 

7 Proof of Homogenization / Gaussian Corrector 
Results 

Here we prove extensions of known results. 
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Figure 10: LDP results for the convolved media. Comparison of empirical rate function 
S{us,i) with theoretical (approximate) rate function Here ^ = 3, \\h\\i = 1, k = 1 and 
e = 1/100 (left) or e = 1/10 (right). Note that the e = 1/100 plot only has i G (0, 0.4)) due 
to the difficulty in sampling the extremely rare event > 0.5. 




Figure 11: LDP results for the convolved media. Comparison of empirical rate function 
£{u^,i) with empirical corrector rate function £{uo + ev,i). Here \\h\\i = 1, k = 1, e = 1/100, 
and ^ = 1 (left) or ^ = 3 (right). In both cases the true rate function is much different than 
that of the corrector. The corrector's rate function performs better when ^ = 3 as can be 
expected by comparison of the pdf for xl and Xi random variables. 

7.1 Proof of theorem 3.1 

Here we prove one-dimensional homogenization results for our media, which has no 
uniform (in lo) ellipticity lower bound and is stationary only in the second variable. 
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7.1.1 convergence 

Solving (1) with = and again with A^ = Aq obtain solutions such as (5). 
Subtracting them we have 

1 1 



Uo{x) - Ue{x) = / F{s) 

Jo 



Ae{s) Ao{s) 



ds 



+ 



{F/Ao) 1 {F/A,) 1 



(iMo) Ao{s) {l/A,} A,{s) 



ds 



= Ii{x)+l2{x). 

We can re-write Ii as 

h{x) = (Fl[o,.]M,) - (Fl[o,.]Mo). 
We are thus motivated to prove the following lemma 
Lemma 7.1. Let H € L°°[0, 1] he deterministic, then 

{{H/A,) - {H/Ao)f < e\\H\\looCA_,. 



(46) 



Proof. Write 
E^{{H/A,)-{H/Ao)f 



(se) 



= £2/ / H{s)H{t)CoY^{se,te,s-t)dtds<e\\H\\looCA-i- 
Jo Jo 

Since < ||/||l2, lemma 7.1 gives us 

E^\\h\\h<s\\f\\l,CA-i, 



A{te, t) 



-E^{l/A} its] 



□ 



(47) 



We now consider I2 . After repeated use of the equality ab—ab= {a — a)b + a{b — b) 
we obtain 



..(.)M(mo)-(FM.))(f ^^^d. 



+ {{l/A,) - (l/Ao)) 



+ 



Jo 



1 



Ao{s) A,{s) 
hiix) + l22{x) + hsix). 



{l/Ao){l/A, 

\ f {F/A,) 



(£/A 

el Jo ^0 



■ds 



Each term is the product of a term similar to Ii and a bounded random variable (a 
priori bounds obtained using the positivity of A). We thus obtain 



^E^||/2||i. < V^E^||/2l||i. + ^E^||/22||i. + ^E^||/23|| 



2 

L2 



The inequality in theorem 3.1 thus follows from (46), (47), and (48). 



(48) 
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7.1.2 a.s. convergence 

Here we use the same decomposition, 

Ue{x) - Uo{x) = Ii{x) + l2l{x) + l22ix) + hsix), 

and show that each term goes to zero a.s. To that end, we notice that every term is the 
product of a bounded (sometimes random) variable, and a term like {H/Af,) — {H/Aq) 
(with H depending on x). It will thus suffice to prove a.s. convergence of this latter 
term. We thus obtain pointwise (in x) a.s. convergence. The a.s. norm convergence 
then follows from an a-priori bound on every realization of — uq and the bounded 
convergence theorem. The a-priori bound follows from (8) (which we assume here is 
independent of e) and (46). In other words, a.s. convergence in theorem 3.1 is a 
corollary of the following lemma. 

Lemma 7.2. Suppose H G L°° is deterministic, then as t ^ oo, we have (almost 
surely P^) 

{H/A^-i) - {H/Ao) 0. 

Proof. The proof is more-or-less a standard trick where we show a.s. convergence 
on a sequence of r values as well as the difference between the sequence values and 
"nearby" values. See section 37.7 in [18]. 

We handle the first convergent term first. Defining 

Y{t): = {H/A,-i)-{H/Ao) 
Xit, s) : = Hit) - {1/A} {t)^ , 

we have 

Now for m G N, a > 0, and m" < r < (m + 1)", 

^y(r) = F(m«) + Z(m^r), Z(m^r) := \ f X{s/T,s)ds. 

Directly from lemma 7.1 we have 

¥.^{¥{111")^] < WHWl^CA-im-", 
and then Chebyshev inequality gives us, for any S > 0, 

P,nK(,„-)i>.i<giiya< ii^iii^g^--'"-' . 
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We choose a > 1 and then 

oo 

^ Pd\Y{m-)\ > 5] < oo. 



So by the Borel-CantelH lemma, y(m") — > a.s. 
As for Z(m",T), we set 



U{m''):= sup \Z{m'',T)\< — \X{s/t, s)\ds, 

m'^<T<(m+l)°- Vm» 



2 /-(m+l)" 

|X(m",r)| < - 

m"<T<(m+l) 

and note that 

\fj\\2 i'(m+l)°' /•(m+1)'" 



\\fj\\2 p(m+ir p(m+ir 

E,|C/(m'^)P < Uf / / C{s-s') 



ds ds' 



< II FTIP r (rn + 1)" - m" 2 ^ r« + l1! 



iL'^A-i ^ II-" lloo'-^A- 

wr m 



Therefore |[/(m")^} ^ m ^, and by Chcbyshcv's inequahty and the Borel-Cantelh 
lemma we obtain U{m"') — >■ a.s. The same conclusion thus holds for Y{t), 
m"- <T<{m + 1)", and therefore for Y{t) = (H/A^-i) - {H/Aq). □ 

7.2 Proof of theorem 3.2 

Here we prove one-dimensional corrector results for our media, which has no uniform 
(in a;) ellipticity lower bound and is stationary only in the second variable. 

7.2.1 One dimensional oscillatory integral 

Here we study the integral 

1 

Ve{x) : = / t)qS) di, qe{t) := -—jr - {l/A} (t), 

Jo ^{t, £ j 

where G{x, s) is deterministic, piecewise continuous in s, and uniformly (in s) Lipschitz 
in X. 

First note that 

/■I r fa~s)/e 

= / G{x,s) / G{x,et + s)Cov^{et + s,s,t)dt 

Jo y-s/e 

Using (9) and dominated convergence, we therefore have 



ds. 



E^ 



^) ^ G{x,tfa^{t)dt, a^it):= r Cov^it,t,q)dq. (49) 

V^J Jo J~oo 



The scaling (in e) and the fact that q^ is mean zero indicate that a central-limit type 
result should show convergence to a Gaussian random variable. This is indeed the case. 
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Lemma 7.3. If G{x,s) is deterministic, piecewise continuous in s, and uniformly (in 
s ) Lipschitz in x, then 



where Wt is a standard Brownian motion. 

The following result allows us to reduce the problem of proving convergence (of a 
stochastic process) to one of studying finite dimensional distributions [5] . 

Proposition 7.1. Suppose (Zn,! < n < oo) are random variables with values in 
the space of continuous functions C([0, 1]). Then converges in distribution to Z^ 
provided that: 

(a) any finite- dimensional joint distribution (Z„(xi), . . . , converges to the 
joint distribution {Z(x,{xi), . . . , Zooixk)) as n ^ oo. 

(b) {Zn) is a tight sequence of random variables. A sufp,cient condition for tightness 
of {Zn) is the following Kolmogorov criterion: There exist positive constants v, 
/3 and 8 such that 



so condition (b-i) is met with v = 2, and (b-ii) is met with (3 = 2 and 5 = 1. 

This means that to prove the theorem we simply need to fix (xi, . . . , and show 




(i) supn>i^ {\Zn{t)\''} < oo, for some t G [0,1] 

(ii) E{\Zn{s)-Zn{t)f}<\t-s\^+', 

uniformly in n> 1 and t,s e [0, 1] . 



Tightness is easily verified. Indeed, (9) and (49) show that 





S-^/^{Ve{x,),... 



dist. 



\v:= {v{xi),. . . ,v{xn)), 




(50) 



Proof. Any finite dimensional distribution has characteristic function 




The above characteristic function may be recast as 




n 



m 
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As a consequence, convergence of the finite dimensional distributions will be proved if 
we can show convergence of 

/,:=/ m{s)^q,{s)ds^Ni:= [ m{s)a{s) dW^, (51) 
Jq V£ Jo 



dist. 



for piecewise continuous m. 

Proceeding, we now show that Ig """""> Ni. We first consider the case of constant 
m=l and Qeis) = q{s, |) = q{^)- To that end, set 



:= / -qi-)ds= / q{s)ds, k = 1, . . . , [e-^ \. 
J{k-\)e ^ ^ Jk-1 



Then /^(.t) = ^j.^^ ^ qk + R{£)i where i?(e) — in (and therefore in probability), 
and hence may be ignored. We therefore consider the limit 

k=i vL^ J k=i 

So we ignore R' and consider the limit of the summation. Following theorem 19.2 in [5], 
we define the sigma fields Tn, T"^ generated by {gfe : A; < n}, {gjfc : A; > n} respectively 
and set 

p„: = sup{|E{r7ir/2}| : m e J^i, r?2 e J^, E {r/^} = 0, E {jyf} = 1, z = 1, 2} . 
Then, so long as Y^=\ Pn <oo, 

VL^ J k=l n=2 

Since p„ < (^(ri — 1), the summability condition on p„ is implied by assumptions 3.1. 
We now show that Indeed, 



E{g2}=E 



1 

q{s) ds 



[ [ Cov^(s-s')dsds'. 
Jo Jo 



Also, using the symmetry of Cov^, we have 

V E {qiqn} = Cov^(s - s') ds ds' = / Cov^(s - s') ds ds'. 

„_o Jo Jl Jo J-oo 



n=2 

Therefore, 



1/2 = / / Cov^(s - s') dsds' = a^. 

Jo J-oo 



This shows (51). 
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To prove the theorem in the case of non-constant m{s), and q = q{s, |) we note 
that if m is replaced by rrih, and q is replaced by qh, giving us 1^^, then 

E^\Ie\x) - = -E^ / / [m{s)q{s, -) - mh{s)qhis, -)] 

£ Jo Jo £ £ (52) 

t t 
X [m{t)q{t, -) - mh{t)qhit, -)] ds dt. 
e £ 

We will choose m/j, (^/j such that the above expectation vanishes in the limit. Split 
[0, 1] into subintervals of size h (with C N) with endpoints tj := jh, j = 1,... , 1/h. 
Now set 

mh{t) := m{tj), and qh{t, ^) := q{tj, ^), for t G [tj,tj+i). 

Taking expectation inside of the integral (52), changing t i-^ et, then t t + s/e we 
have (similar to (49)) 

1 /•oo 



pi poo 

E^|//-7,|2= / / Q,^h{s,t)dtds, 

Jo J-oo 



where \Qe,his,t)\ < Cr{t) (with T from (9)). Using the continuity (on the diagonal) of 
Cov^, and the piecewise continuity of G, we also have Qe,h — >■ as (e, /i) (0,0) (in 
any way whatsoever, for a.e. (s,t)). Therefore, it suffices to prove the lemma 7.3 with 
m/i, Qh replacing m, q. 

We proceed to split the integral defining up. 

Ie -= rnh{t)qh{t, -) dt = / m{tj)q{tj, -) dt. 

Jo £ ^.^0 J[tj,tj+i) £ 

Each subintegral is handled exactly as before, yielding 

m{tj) f q{tj, -)dt^ j mitMtj) dWf 

It remains then to show that the limiting Gaussians (above) are independent. We 
do this in the case where = 1^,1 + 1^,2 is split into two intervals, the general case 
following by induction. To that end, we show that for all fci, A;2 G M 



f : = |E; 

Define 



|gi(fcl/l,e+fc2/2,e)| _ je'^^l^l.^ | E^ |e^*^2/2,.|| ^ Q_ (53) 



Ix+T] 

Then the first term in £ decomposes as 



Ph ■ = n 4=Ab QI2 ■■= Ie,2 - Pi,. 

Jx+Ti y£ ^ 
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The first term in the decomposition is small since 
As for the second term, 

- [e^ { (^e^'^i^M - 1) e^^^i^M } E^ |e**^2/.,2 1" _ 

The first bracketed term is < 99(77/6) by our mixing condition (3.1), and the second is 
< 7^1/2 ^ manner similar to (54). Therefore, 



as £ ^ 



for say rj = e^/^. 

We have thus shown 



^ [' mh{t)ah{t)dWt, 
Jo 



J dist 



and this completes the proof. □ 

We also prove a complementary lemma that allows us to deal with the lack of a 
uniform (in oj) ellipticity lower bound. 

Lemma 7.4. Suppose A satisfies the hypothesis of theorem 3.2. Then with 

s 1 1 

q{s,-) ■ = 



we have for any H G -L°°([0, 1]), 

f := E^ I (^j^^ H{s)q{s, J) | < Q j^^ H{s)H{t)p^'\^^) ds dt^ 



Proof. 

£< I \H{si) ■ ■ ■ H{se)\ IEj \q{si, ^) • • • qis^, ^)|| ds. 

The expectation can be broken up using Cauchy-Schwartz into a product that looks 
like 

(E5 {Xfxl} E^ {Xlxi} E^ {Xlxl})'/' . 
For each term [Xf Xf^-^^ , the mixing condition, and our bound on 6*'^ moments 



gives 



E^|,(.,J)3g(t,^)3}</,(^)Q. 



The result follows. □ 
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7.2.2 Asymptotic expansion of the solution 

We now complete the proof of theorem 3.2. Starting from (6), (7) we write 

Ue{x)-Uoix) 1 . X . 1 „ / X 

Vi = ^-e(x) + ^i?.(x). 

The convergence of £~^^^Ve is assured by lemma 7.3. We now show that £~^/'^Re is tight 
and converges pointwise to zero. Then, proposition 7.1 shows that £~^/^i?£ '^^^*'> in 
the space of continuous paths. 

To that end we write £~^/'^R{x) as a sum of terms of the form 



Jo 



1 



(. 



1 



1 



10 \Ms) Ms] 

for a constant random variable Z, with either 

\Z\ < C'^iXff, \Z\ < C'^\XfYf\, \Z\ < C'^\Xf\{l/A,), or \Z\ < C'^Yf. 

We first show that [{e~^/'^{Re{y) — Reix)))'^} < e\y — xp, meeting condition (b-ii) 
of proposition 7.1 with /? = 3, (5 = 1. Then choosing x = 0, we have (since /^^(O) = 0) 
£~^/'^Rs[y) -> in L^{Q,,'P^). This meets condition (b-i) of proposition 7.1 with = 3 
and also shows that all finite dimensional distributions converge to the zero vector as 
well (meeting condition (a) the proposition). 
Using Holder, with \/p + \/p' = 1, and r > 



y 1 



X Be{x) 



rp 



We choose r = 2, p = 3/2, p' = 3 to get 



Now due to lemma 7.4, 
IE. 



Be{x) 



ds 



1/3 



ry 1 
Jx 'bJx) 



ds 



1/3 



^ W llPllool 



— X\ 



Also, using Holder and lemma 7.4, E^ {l-Z^P} is < £ in all cases of Z. We therefore 
have 



E. 



RM-Re{x) 



< eCl\y-xf. 



Acknowledgment 

The authors would like to thank Mark Adams and Yu Gu for many helpful discussions. 
This work was supported in part by NSF grant 0904746 and NSF RTG grant DMS- 
0602235. Any opinions, findings, and conclusions or recommendations expressed in 
this material are those of the author(s) and do not necessarily reflect the views of the 
National Science Foundation. 



38 



References 



[1] G. Bal. Central limits and liomogcnization in random media. Multiscale Model. 
Simul, 7(2):677-702, 2008. 

[2] G. Bal, J. Garnier, S. Motsch, and V. Perrier. Random integrals and correctors 
in homogenization. Asymptot. Anal, 59:1-26, 2008. 

[3] G. Bal and W. Jing. Corrector theory for msfem and hmm in random media. 
submitted, 2010. 

[4] G. Bal and W. Jing. Homogenization and corrector theory for linear transport in 
random media. To appear: Disc. Cont. Dyn. Syst. B, 2010. 

[5] P. Billingsley. Convergence of Probability Measures. John Wiley and Sons, New 
York, 1999. 

[6] Bousquet O. Bouchcron, S. and G. Lugosi. Concentration inequalities. Advanced 

lectures in machine learning, pages 208-240, 2004. 

[7] A. Bourgeat and A. Piatnitski. Estimates in probability of the residual between the 
random and the homogenized solutions of one-dimensional second-order operator. 
Asympt. Anal, 21:303-315, 1999. 

[8] A. Bourgeat and A. Piatnitski. Approximations of effective coefficients in stochas- 
tic homogenization. Ann. I. H. Poincare, 40:153-165, 2004. 

[9] W. Bryc. On large deviations for uniformly strong mixing sequences. Stochastic 
processes and their applications, 41:191-202, 1992. 

[10] J. Bucklew. Introduction to rare event simulation. Springer, 2004. 

[11] R.M. Burton and H. Dehling. Large deviations for some weakly dependent random 
variables. Statistics & probability letters, 9:397-401, 1990. 

[12] R. E. Caflisch. Monte carlo and quasi-monte carlo methods. Acta Numerica, pages 
1-49, 1998. 

[13] A. Dembo and O. Zeitouni. Large deviations techniques and applications. Appli- 
cations of mathematics. Springer, 1998. 

[14] A.B. Dieker and M. Mandjes. On asymptotically efficient simulation of large 
deviation probabilities. Adv. Appl. Prob., 37:539-552, 2005. 

[15] G. Folland. Real analysis, modern techniques and their applications. John Wiley 
& Sons, 1984. 

[16] R. G. Ghanem and P. D. Spanos. Stochastic finite elements: a spectral approach. 
Springer- Verlag, New York, 1991. 

[17] V. V. Jikov, S. M. Kozlov, and O. A. Oleinik. Homogenization of differential 
operators and integral functionals. Springer- Verlag, New York, 1994. 

[18] M. Loeve. Probability Theory II. Springer, fourth edition, 1978. 

[19] Owhadi H. Lucas, L.J. and M. Ortiz. Rigorous verification, validation, uncer- 
tainty quantification and certification through concentration-of-measure inequali- 
ties. Comput. Methods Appl. Mech. Engrg, 197:4591-4609, 2008. 

[20] G. Pavliotis and A Stuart. Multiscale Methods: Averaging and Homogenization. 
Springer, 2008. 

[21] C. Robert and G. Casella. Monte Carlo Statistical Methods. Springer, 2004. 



39 



